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И ОБРАБОТКА РЕЗУЛЬТАТОВ ПРИ ПОМОЩИ 
ПРОГРАММНОГО СРЕДСТВА ЕО$УРРОВТ 


Описываются возможности программного средства ЕОЗиррогЕ разработанного ав- 
тором и предназначенного для поддержки принятия решений в задачах, требую- 
щих планирования эксперимента и статистической обработки результатов. Приво- 
дятся научные результаты некоторых практических применений программы. 
Ключевые слова: планирование эксперимента, дисперсионный анализ, описание 
поверхности отклика, оптимизация процесса. 


Введение. Существует множество программных средств, предназначенных 
для автоматизации вычислений и осуществляющих поддержку принятия 
решений в разнообразных задачах, связанных с планированием экспери- 
мента. Однако практически все они являются разработками зарубежных 
университетов или фирм, не поддерживают русскоязычный интерфейс и 
являются платными программными продуктами [4]. Поэтому является акту- 
альной задача создания собственной программы для поддержки планиро- 
вания экспериментов. 

Реализованные методы. Созданная ранее автором функционирующая 
версия программы ЕОбирро“ (Ехрейтета! Ор{ттхаНоп Зиррог, Система 
поддержки экспериментальной оптимизации), позволяет строить линейные 
и квадратичные регрессионные модели (содержащие до 10 факторов) ис- 
следуемых процессов, а также проводить статистическую обработку полу- 
чаемых результатов. 

Для построения линейной модели пользователю предоставляется 
возможность выбора между полным и дробным факторным экспериментом 
(ПФЭ и ДФЭ соответственно). В случае выбора полного факторного экспе- 
римента в линейную модель можно дополнительно включить произвольные 
эффекты взаимодействия факторов. При проведении факторного экспери- 
мента используется схема с одинаковым количеством повторных опытов. 
Эти дополнительные опыты позволяют оценить дисперсию единичного из- 


ы 
мерения ее : 
№ т у М ь т 
[У„ ^ У. 5 У» 
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где 1= \..№М - номер строки матрицы планирования; и-= \..т - номер по- 
вторного опыта; 5, - дисперсия отклонения повторных опытов в 
каждой точке плана от среднего значения отклика в этой точке; 
У - матрица № т откликов; число степеней свободы ое равно 


{= Мт- ). 
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Коэффициенты регрессионной модели рассчитываются также на 
основании средних значений У,, что обеспечивает высокую точность их 
оценки: 


т . ЕЕ ТОЙ (2) 


где [ - количество коэффициентов в модели. 
Статистический анализ линейной модели состоит в следующем. 


м у 
Сначала проверяется гипотеза об однородности дисперсий 5, по критерию 


Кохрена. Затем при помощи [’-критерия Фишера проверяется гипотеза об 
адекватности модели и, наконец, для каждого рассчитанного коэффициен- 


та В ‚ можно проверить его значимость по # -критерию Стьюдента [1]. Дан- 


ный набор критериев позволяет достаточно полно оценить надежность по- 
лученных результатов эксперимента. 

Пример 1. Использование ПФЭ для исследования взаимодействий 
факторов в «игрушечном гольфе». Рассмотрим эксперимент по забрасыва- 
нию мяча на заданное расстояние при помощи специальной установки, на- 


зываемой «игрушечный гольф» [3]. Управляемые факторы: № - длина 
биты (или клюшки); Х. - угол отклонения биты перед ударом; »Х. - масса 


дополнительного груза биты. В качестве отклика выступает расстояние У, 
которое преодолеет мяч после удара. Необходимо исследовать взаимодей- 
ствие между факторами и построить регрессионную модель для предсказа- 
ния результатов произвольного удара. 
Методы испытаний. Для решения поставленной задачи подходит полный 
факторный эксперимент, позволяющий получить неполную квадратичную 
модель, содержащую линейные члены и эффекты парных взаимодействий. 
Интервалы варьирования факторов и некоторая дополнительная информа- 
ция представлены в табл.1. Сгенерированный ЕОбиррой план эксперимен- 
та и результаты непосредственных измерений представлены в табл.2. 
Таблица 1 
Общие данные об эксперименте 


т «А С 


| Тип плана ____| плана Полный факторный варьирование факторов на 2 уровнях 
Модель Е 
еб еовЕ еЕФа 
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Таблица 2 
Матрица планирования и результаты непосредственных измерений 


ПС Е О Е Е О Е Е БТ С 
НЕ ТЯ ЖЕ Е Е О ТН ЕН НЕВЕ ПИРА 





Результаты эксперимента и их обсуждение. В результате опытов 
были получены коэффициенты уравнения регрессии. Проведен статисти- 
ческий анализ модели в целом и ее коэффициентов в отдельности. Ре- 
зультаты сведены в табл.З3. 
Таблица 3 
Результаты статистического анализа эксперимента 


Коэффи-ци- | Степень сво- Оценка Ё -экспер. Ё -крит. ( Вывод 
енты боды 
РВ =0.05) 
\ 
ы 
у 


значим 
у 


10.58 
| 


Модель Е -эксп. Е -крит. ЕЕ 
Р =0.01) 
кватна 





Получена следующая регрессионная модель в кодированных едини- 
цах: 
у= ТУ.Л+ \1.Ах + ох, Л.А Ах, А Ихх, + 
Е ИЖ, 


122 


Раздел «Управление, вычислительная техника и информатика» 








В натуральных единицах: 
у= 1.8 - МУХ - 4 Тх,- 


АЕ Е ТЕ. 
из которой видно, что из всех парных эффектов взаимодействия наиболее 
сильным является взаимодействие хх, (длина биты и масса груза). Все 


вычисления проводились два раза: первый — при помощи авторской про- 
граммы ЕОбиррой и второй — при помощи коммерческого продукта 
Рез!ап-Ехрек {па|-версии 7.1.4 фирмы ЗаЕЕазе, 1пс. Результаты полностью 
совпали. 

Описание экстремальной области. После получения линейной модели, 
описывающей локальную область факторного пространства, исследовате- 
лю надлежит принять решение о дальнейшем ходе исследования; следует 
ли оптимизировать процесс, двигаясь в направлении градиента (траекто- 
рию автоматически рассчитает ЕОбиррог), либо область экстремума уже 
достигнута и необходимо построение более точной квадратичной модели. 
Чтобы выяснить это в ЕОбиррог!, существует специальная функция оценки 
кривизны поверхности отклика в области эксперимента. Для этого пользо- 
вателю предлагается поставить несколько дополнительных опытов в цен- 
тре плана. Мерой кривизны поверхности отклика может служить разность 
среднего значения этих откликов и свободного члена линейной модели [2]. 

Построение квадратичных моделей также легко осуществимо при 
поддержке ЕОбиррок. Для этого в текущей версии имеется поддержка по- 
строения ортогональных центральных композиционных планов (ОЦКП), в 
которых факторы варьируются на пяти различных уровнях. Достаточно 
только указать количество факторов, количество повторных экспериментов 
и количество дополнительных экспериментов в центре плана, и система 
автоматически рассчитает величину «звездного плеча» @ , а также смеще- 
ние для квадратичных членов. Также можно достроить существующий план 
ПФЭ или ДФЭ до ОЦКП. 

Обычно на практике проведение экспериментов является дорого- 
стоящей и трудоемкой операцией, поэтому при планировании всегда стре- 
мятся уменьшить количество точек плана. Чтобы достичь этого без суще- 
ственных потерь точности, применяется специальная схема, согласно кото- 
рой повторные опыты, необходимые для оценки дисперсии единичного из- 
мерения, ставятся только в центре плана [5]. Использование данной схемы 
иллюстрирует следующий пример. 

Пример 2. Лрименение ОЦК! для оптимизации процесса очистки 
воды. Вода, полученная после процесса осушения баков для хранения бен- 
зина, подвергалась очистке в пробной установке для исследования эффек- 
тивности совместного влияния коагулянта (сульфата алюминия) и катионо- 
вого флокулянта (полиэлектролита) группой ученых из Национального тех- 
нического университета в Афинах [6]. 
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Методы испытаний. Для построения квадратичных регрессионных моде- 
лей трех параметров оптимизации (мутность, наличие твердой взвеси, на- 
личие нефтепродуктов) использовалось композиционное планирование. В 
качестве факторов были выбраны концентрации коагулянта, флокулянта и 
уровень рН. Автор настоящей статьи использовал полученные исследовате- 
лями экспериментальные данные при апробации программы ЕОбиррой для 
проверки вычисленных статистических оценок. Пример ограничивается 
лишь одним параметром оптимизации — уровнем замутненности воды, под- 
вергшейся очистке. Общие данные представлены в табл.4 и 5. 

Таблица 4 


Общие данные об эксперименте 


| Фактор | Название | Единица | -1.414 |1 [0 | +1 | аам | 


\ 
у 
ЕЯ Зи БАБ АН 


Отклик _| Название | Единица | | 


и Е 


ОЦКП (варьирование факторов на 5 уровнях 
> у у \ 
У= Ь. + Вх, + Бх, + Вх. + Б,хх, + Быхх, + Быхх. + Вх + Вх, + Вых, 
Количество опытов в центре плана 





Таблица 5 
Матрица планирования и результаты непосредственных измерений 


ВЕТ ВЕ ЧИ ЕН ВЕСИ 
ра | 69| мо | 141 | 
па [евро |141 | 
кр [256 [во о | 
[689 | рао|ос 
ЕВ о 
О Е ТЕ ВЯ ВЕ ЗИЕ ОЕЗНЕ 
риа ос 
С Е Е ОС О ЕЙ 


Результаты эксперимента и их обсуждение. После проведения опы- 
тов были получены коэффициенты квадратичного уравнения регрессии. 
Проведен статистический анализ модели. Результаты представлены в 
табл.е6. 
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Таблица 6 


Результаты эксперимента 


Г ВЫ [В [вы | 


26.95 
ЕВЕ СЕНЕ 


ДЫ 


у 5 267.45 
5 д 
у 3 21.41 
5 вост 
“крит. (Р =0.05) 
нии 


12.49 9. 


Особый интерес представляют графики полученной поверхности 
отклика, представленные на рисунке, определяемой уравнением (в кодиро- 
ванных значениях факторов): 


ПЕ АА АТО о: 
+4 Ух. Ах + УЛ + 4. 


На каждом из трех графиков отсутствующий фактор зафиксирован на 
основном уровне (0.0). 








Зависимость мутности воды (У в %) от дозы 


полиэлектролита (Х. ), сульфата алюминия ( 


Х, ) и уровня РН (.Х..) в кодированных зна- 





чениях 
Из рисунка хорошо видно, что мутность воды в значительной 


степени зависит от дозы полиэлектролита и в меньшей степени от уровня 
рН и дозы сульфата алюминия. Кроме того, мутность имеет минимум около 
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точки (-. м ме .У) на уровне около 25%. Результаты, полученные при 


помощи программы ЕОбиррром, хорошо согласуются с результатами 
исследователей, изложенными в [6]. 

Выводы. Программное средство ЕОбиррой планируется в будущем расши- 
рить дополнительными возможностями, такими как поддержка ротатабель- 
ных планов, планов Бокса-Бенкена, функциями разбиения плана на блоки, 
вычисления р-значений для проверки гипотез и др. Наряду с улучшением 
средств визуализации полученных результатов всё это позволит суще- 
ственно повысить эффективность применения ЕОбирррой. 
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